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Abstract. In the present paper we present JHK photometric results of the young embedded cluster NGC2316. We construct 
the cluster radial profile from which we determine a radius of 0.63 pc. We find 189 ±29 cluster members in an extinction 
limited sub-sample of the survey, 22 ± 19 of which are possibly substellar. An average extinction of 4.5 visual magnitudes is 
derived using (H - K) colours of control fields. This extinction is due to the presence of residual parental molecular cloud. 
NGC 2316 presents 16% source fraction of excess emission which is consistent with other results from clusters with an age of 2 
- 3 Myr This age is consistent with the distribution of sources in the colour-magnitude diagram when compared to theoretical 
isochrones, and the overall shape of the cluster KLF The substellar population of the cluster is similar or smaller than that 
observed for other embedded clusters and the stellar objects dominate the cluster membership. 
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1. Introduction 

Embedded clusters may be the fundamental units of star forma- 
tion (Lada & Lada, 2003). They contain statistically significant 
samples of young stars of similar age and composition span- 
ning a wide range of mass, providing excellent laboratories for 
investigating a number of important issues in star formation 
such as the form and universality of the initial mass function 
(IMF) and the frequency and lifetimes of protoplanetary disks. 
Surveys for near infrared (NIR) excess in embedded clusters 
can provide knowledge of the disk frequency and evolution 
because a protoplanetary disk is easier to detect than a plan- 
etary system with a similar mass of solid material (Beckwith 
& Sargent, 1996). In addition, embedded clusters provide the 
smallest spatial scale for investigating the nature of the IMF, 
its form and universality; their mass function is in fact an ini- 
tial mass function since young clusters have not lost significant 
numbers of members to either dynamical or stellar evolution. 
NGC 23 16 is a young partially embedded cluster with coor- 
dinates (a,5)(J2000) = (6^59'"40%-7°46'36") at a distance 
of 1.1 kpc according to FelH et al. (1992) and Hodapp (1994). 
There is an HII region centered on a B3 ZAMS star (Noguchi et 
al., 1993; Fukui et al., 1993), associated with the IRAS 06572- 



0742 source. The UV radiation is apparently creating a dense 
photodissociated region (Ryder et al., 1998), showing some 
physical signatures of ongoing star formation, such as H2O 
masers (Felli et al., 1992) and CO outflow with spectral index 
characteristic of optically thin free-free emission (Beltran et 
al., 2001). NGC 23 16 was first uncovered by Parsamian (1965), 
followed by other studies (e.g. Hodapp, 1994), however certain 
parameters are still not well established. To better characterize 
this cluster we obtained deep multi wavelength NIR images of 
the region to better constrain the cluster's age and spatial dis- 
tribution, as well as provide estimates of membership number, 
fraction of NIR excess sources and brown dwarf population. 

2. Observations 



Fig. 1. (See attached file Figl.jpg) NGC 2316 JHKs combined 
grey scale image from the NTT data.The scale is stretched to 
enhance the low level emission due to residual cloud material. 
The large majority of sources in this image are not visible at 
optical wavelengths. The field size is 5'x5'. 



Send offprint requests to: P. S. Teixeira 
* Based on observations carried out at ESO, La Silla, Chile. 



The observations were made on the 8th and 9th of March 
1999 using the Sofl infrared camera (/(1.25//m), //(1.65jum), 
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and /Ti (2.162yL/m) bands) on the 3.5 m NTT telescope in La 
Silla, Chile. We took 45 images per band, with a total integra- 
tion time of 900 s. Two control fields were taken also in JHKs 
bands, one located 20' N of the NGC2316's center, the other 
located 13' S, with a total integration time of 90s. The field of 
view was aproximately 5' x 5', with a pixel scale of 0.29"per 
pixel. The average FWHM of point sources in the final image is 
about 0.54". In Fig. 1 we can observe that NGC 23 16 presents 
a residual cloud material, nicely visible as a bubble centered 
on the cluster's core. This faint nebulae, surrounding the B 
star in the center of the cluster, is indicative of the early age 
of this cluster, less than 5 Myr when compared to NGC 2362 
( |Moitinho et al., 2001t which shows no nebulosity. 

3. Data analysis 

The data reduction for NGC 23 16 cluster is done using the 
IRAF' DIMSUM (Deep Infrared Mosaicing Software) pa- 
ckage, to reconstruct dithered exposures taken for J, H and 
Ks filters. Astrometry is conducted with ESO's (European 
Southern Observatory) SKYCAT tool and the IRAF routines 
IMTRANSPOSE, CCMAP and CCS ETWCS. Infrar ed sources 
are identified using SExtractor (Bertin & Arnouts, 1996.) and 
checked through visual inspection. A FWHM of 0.5"is used, 
along with a gaussian filter of 3 sigma, a detection threshold of 
3 sigma and a deblending threshold parameter of 64 with a mi- 
nimum contrast of 0.0005 counts. We use the SExtractor output 
to remove approximately 80 galaxies from the cluster and con- 
trol fields K-hwA data. Photometry is done with the routine 
PHOT of the IRAF's PHOTCAL package, using an aperture 
radius of 5 pixels (0.73") for the cluster and 6 pixels (0.87") 
for the two control fields with appropriate aperture corrections 
per band. Finally, for the photometric calibration we reduced 9 
exposures of the standard star 9118 from the Persson catalog 
l lPersson et al., 1998j , determining the zero point by using the 
filter transformation equations in the Sofl's camera web page^. 
We converted our Ks photometric results to K fotometry using 
these equations to better compare with other similar studies on 
young clusters. Comparing our photometry with 2MASS pho- 
tometry^, we find an average difference in 7, H and K magni- 
tudes of 0.04, 0.02 and 0.04, respectively and an average dif- 
ference in {H - K) colours of -0.01 and in { J - H) of 0.03, so 
no further calibration is needed. 

4. Results 

In the present JHK photometry survey of NGC 23 16, we de- 
tected 1067 sources at 7-band, 1229 sources at //-band , and 

' Image Reduction and Analysis Facility (IRAF) is distributed by 
NOAO, which is operated by AURA, Inc., under contract to the NSF. 

^ http://www.ls.eso.org/lasilla/Telescopes/NEWNTT/sofi/. 

^ This publication makes use of data products from the Two 
Micron All Sky Survey, which is a joint project of the University 
of Massachusetts and the Infrared Processing and Analysis 
Center/California Institute of Tecnology, funded by the National 
Aeronautics and Space Administration and the National Science 
Foundation. 



Fig. 2. (See attached file Fig2.jpg) /sT-band radial profile of 
NGC 23 16. The error bars represent the V/V statistical error in 
each bin. The dashed line coiTesponds to the background den- 
sity determined from averaging the reddened ((Av)=4.5 mag- 
nitudes) control field. Two fits are plotted, one using the King 
(1962) model, the other an inverse radius profile. From this we 
estimate the cluster radius to be 120", or 0.63 pc. 

Fig. 3. (See attached file Fig3.jpg)(// - K) vs. ( J - H) colour- 
colour diagram for NGC 2316 sources. The giant branch and 
main sequence from Bessell & Brett (1988) are shown, along 
with the Classical T-Tauri locus from Meyer et al. (1977). The 
reddening vector was determined using the values of extinction 
from Rieke & Lebofsky (1985). Photometric errors are illus- 
trated at the upper right. 

1251 sources at K-hwd. The completeness limits of the ob- 
servations are 18.0, 18.5 and 20.0 for the K, H, and 7-bands, 
respectively, in the cluster field, and 17.5, 17.5, 18.0 for an a- 
veraged control field. These are estimated from the luminosity 
functions, where we consider our sample complete up to the 
bin imediately before the maximum. For the data analysis the 
control field completeness limits were used and no extrapola- 
tion was made to account for sensitivity diferences between the 
two fields observed. 

4.1. Radial profile 

Figure 2 presents the K radial density profile, where only the 
sources brighter than the completeness limit were included. 
The profile is centered on the peak of the stellar density distri- 
bution. NGC 23 16 presents a regular distribution falling steeply 
and merging with the extincted background at 120"(0.63pc). 
The stellar surface density distribution exhibits more of a cen- 
trally condensed and relaxed structure than a hierarchical type 
structure (e.g., Lada & Lada, 2003). Within the 120" radius, the 
number of cluster members is 158±20. Assuming a spherical 
geometry and an average mass per star of 0.3 M©, we determine 
the stellar volume density to be 45 M© pc"^. Two models are fit- 
ted to our radial profile as shown in the figure. The King model 
(King, 1962 1 is given by the expression n(r) - /o/[l -i- ('"/'"c)^], 
where fo is the core concentration at zero radius and is 
the core radius. This function agrees sufficiently with our data 
(X^ = 0.81) to provide an estimate of the core radius and den- 
sity. The parameters are determined to be 300 stars arcmin"^ 
(72 Mopc"-') and 10" (0.05 pc) respectively for fo and r^. An 
inverse radius model (n(r) - -5.77 + I r, - 0.94) 

also describes fairly well the density distribution. We also note 
some oscilations in the density distribution that can be of sta- 
tistical origin. 

4.2. Colour-colour diagram 

The infrared JHK colour-colour diagram of NGC 23 16 is plot- 
ted on Fig. 3, considering only sources brighter than the com- 
pleteness Umits for each band and located inside the cluster 
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Fig. 4. (See attached file Fig4.jpg)/sr vs. {H - K) colour- 
magnitude diagram for the entire survey. The data are com- 
pared with the 1 and 5 Myr isochrones of Baraffe et al. (1998) 
at a distance of 1.1 Kpc. A sample is selected for sources with 
Ay < 10; the reddening slope is determined from Rieke & 
Lebofsky (1985). 



radius, as determined in the previous section. There is a con- 
siderable spread of stars along the reddening band, which cor- 
roborates with the previous statement that the cluster is par- 
tially embedded. {J - H) and (H — K) colours are also calcu- 
lated for the two control fields for comparison. The latter fields 
show almost no reddening. Subtracting this contamination we 
find 16%± 3% of cluster members with NIR excess emission 
characteristic of young stars with circumstellar disks (Lada & 
Adams, 1992). 

4.3. Colour-magnitude diagram 

The infrared colour-magnitude diagram for NGC 2316 is pre- 
sented in Fig. 4 for sources with magnitudes above the com- 
pleteness limits. We compare source locations with the theo- 
retical isochrones of 1 and 5 Myrs from the Baraffe et al. 
(1998) non-grey evolutionary models at a distance of 1.1 Kpc 
dPelh et al., 1992) . There is some spread of sources to the left 
of the isochrones due to foreground contamination. From the 
diagram, we observe a much wider spread of sources to the 
right of the isochrones, a result of extinction produced by the 
associated molecular cloud. In order to derive the cluster KLF, 
we construct a statistically more complete mass and extinction 
sample by including only objects with Ay < 10 magnitudes 
and mass > 0.04 Mq, for which we are essentially complete. 
This diagram is not a good discriminant for the cluster age so 
we assume an intermediate age of 3 Myrs for the construction 
of our mass limited and extinction limited sample. Within the 
errors this sample is also valid for cluster ages of l-5Myr. 

4.4. Luminosity function 

The absolute K luminosity function of the cluster (KLF), ad- 
justed for extinction, is plotted in Fig. 5, where the hydro- 
gen burning limit is indicated by a vertical line. Here we 
only use sources from the sample determined from the colour- 
magnitude diagram. This KLF is determined by subtracting 
a reddened control field from the original reddened on field. 
The off" fields are averaged into a final control field and scaled 
so that the control field and the on field have both the same 
area. The extinction is determined from the average colour ex- 
cess of the NGC 2316 field as compared to the control field 

(A(H - K) ^ {(H - /:)cluster) - {(H - /:)contromeld)) USing 

sources with no NIR excess emission. This average colour ex- 
cess is used to determine the extinction following the Rieke 
& Lebofsky (1985) law and is found to be (Ay) = 4.5 magni- 
tudes with a standard deviation of approximately 5 magnitudes. 
We compare this luminosity function with those of Trapezium 
dMuench et aLr2002,> and IC 348 ( .Lada & Lada, 1995^ , where 



Fig. 5. (See attached file Fig5.jpg)Dereddened KLF of 
NGC23 16, compared to the KLFs of the Trapezium and IC 348 
clusters. The vertical dotted line indicates the hydrogen burn- 
ing limit (HBL), jBarafi^e et al., 1998> for NGC 2316 assuming 
a distance of 1 . 1 kpc and an age of 3 Myr. 

we find a striking similarity between the shapes of NGC 2316 
and IC 348 KLFs. We estimate a membership of 189 ± 29 and 
a number of objects beyond the HBL of 22 + 19. This is a very 
crude estimate because although we are sensitive to masses as 
low as 20 Mjup, the background contamination dominates the 
statistical errors for the substellar KLF bins at NGC2316's dis- 
tance. 

5. Discussion and summary 

Circumstellar disks emit at NIR wavelengths. Studying the NIR 
excess emission of young cluster sources provides thus infor- 
mation on the frequency of circumstellar disks. Our estimate 
of the NIR excess emission for NGC 23 16 according to the 
colour-colour diagram is 16% +3%. 

We conclude that NGC2316's members span a wide range of 
masses, from B-stars down to the substellar regime (40Mjup), 
and is significantly embedded ((Ay) = 4.5) in its parental 
molecular cloud. We estimate a size of 0.63 pc for NGC 2316 
from the radial profile assuming a distance of 1 . 1 kpc. Using 
an extinction limited sub-sample of the survey, we derive a 
membership for the entire cluster of 189 + 29 objects, of which 
22 + 19 are beyond the HBL. 

Regarding the crude estimate of the substellar population we 
note that studies of clusters at distances around and beyond 
1 kpc will be hampered by an overwhelming background con- 
tamination at the faintest magnitudes (that cannot be reduced 
by deeper observations). Nevertheless, the estimated fraction of 
substellar objects in NGC 2316 (2-22 %) is comparable, within 
the errors, to that derived for the Trapezium cluster (20-25 % 
Muench at al., 2002) and IC 348 (14-20 % Muench et al., 2003) 
from a similar KLF analysis. This is additional evidence that 
brown dwarfs do not make a significant contribution to the 
overall cluster, either by numbers or mass. In other words, the 
population of these clusters is dominated by stellar rather than 
substellar members. 

The NIR excess emission obtained for NGC 2316 from JHK 
photometry is likely an underestimate. Significant emission 
of a circumstellar disk is in the L-band (3.4 /vm), arising 
from warm circumstellar material within a few stellar radii, 
so JHKL photometry would provide a more complete cen- 
sus of excesses, ULiuetaL, 2003j . IC 348, a cluster of appro- 
ximately 3 Myr, presents a fraction of NIR excess sources of 
21% using JHK photometry and 65% using JHKL photome- 
try, dHaisch et al., 2001 Lada et al, 2000.1 . On the other hand, 
NGC 2316 has a bright UV source, a central B3 star, which 
might accelerate the dissipation of circumstellar disks by pho- 
toevaporation and consequently decrease the number of ob- 
served sources with NIR excesses. 

The comparison of the NGC 23 16, IC 348 and Trapezium KLFs 
in section 4.4 shows good agreement between their shapes; in 
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particular, the first two are very similar (apart from the num- 
ber of sources), as they roughly show the same broad peak 
at the same absolute magnitudes. Assuming an universal IMF, 
(e.g. lLada & Lada, 2003) , we infer that NGC 2316 has approx- 
imately the same age as IC 348, about 2-3 Myr (Lada & Lada, 
1995; Muench et al., 2003), but is more evolved and older than 
the Trapezium. 
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